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Direct Distance Measurements to SN 2009ip 
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ABSTRACT 

We show applicability of our new method, Dense Shell Method (DSM), for measuring 
astronomical distances using an example of supernova SN 2009ip, which had several 
outbursts. The method is based on the original Baade idea. Distance to this supernova 
is well measured via Cosmic Distance Ladder and generally assumed to be 20.4 Mpc. 
Our method directly gives a very close result, namely we find that the median distance 
is rs 20.1 ± 0.8 (68% CL) Mpc to SN 2009ip assuming the most reasonable values for 
its parameters. 
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1 INTRODUCTION 

We l|Blinnikov et alj|2012h have introduced a new method 
for measuring cosmological distances, which we suggest to 
call the Dense Shell Method (DSM). The method relies on 
the formation of an expanding dense shell in SN Iln and 
allows to find a linear size of a supernova shell in absolute 
units and distance to it without h elp of the cosmo l ogica l 
distance ladder. Using this method iBlinnikov et al.l (|2012l ) 
have calculated the distance to SN 2006gy. But the result 
contains quite a large error due to uncertainty of the in- 
terstellar extinction and the color temperature. Moreover, 
there is a small number of observations of the rising part of 
the light curve when the shell is believed not to fragment. 
In this paper, we examine applicability of our method for 
the case of the supernova SN 2009ip, which exploded again 
in 2012. We show that the accuracy of the DSM distance 
determination in this case is much higher. 



2 DIRECT DISTANCE DETERMINATION BY 
DENSE SHELL METHOD 

Supernovae of type Iln are obse r ved at very large r edshifts 
!|Cookell2008l ; ICooke et ai]|2009l ; iMoriva et al.|[2012l ). They 
do not ente r the coasting free expan sion phase, t herefore nei- 
ther EPM (jKirshner fc Kwanlll974l ) nor SEAM l|Baron et all 
|2004 ) techniques are directly applicable. Luminous Blue 
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Variables can eject surface layers several times, filling cir- 
cumstellar space by an expanding cloud of relatively dense 
matter. The kinetic energy of those ejecta may be sub- 
stantially lower than t hat of the genuine supernova (see 
iHeeer fc Wooslevl l2002h . Narrow lines in the spectra of 
SN Iln imply that the matter of first ejections has velocity 
an order of mag nitude lower than in the case of the ordinary 
supernovae (see iGrasberg fc Nadezhin|[l98r3) . 

SN 2009ip has produced several outburst which were 
observed directly. The physical mechanisms of t hose ejec- 
tions can be different jHeger fc Wooslevl |2002| ; IChevalierl 
l2012i ; ISoker fc Kashill2012i ). The time between the outbursts 
can greatly vary. If this time is large enough, a "dip" can 
be formed between the star and the locus of the maximum 
density of the cloud. If one imagines that the LBV giant sur- 
rounded by such a cloud from the primary explosion would 
outburst again as a genuine supernova or produce a simi- 
lar but high-velocity outburst the new ejected material will 
have much higher velocity. So in the first days after the 
shock breakout, ejecta would fly through the surrounding 
low-density "dip" of the cold and thus transparent cloud. 
The spectra would contain only the broad lines born in 
the high-velocity ejecta. Ejecta would propagate and col- 
lide with the high-density cloud layers. Shock waves would 
run to both sides of the locus of the collision, forward (into 
the cloud) and backward (into the ejecta). 

In what follows we consider only spherically symmetric 
outbursts. Our calculations show that the velocity of these 
shocks relative to ejecta is low. This is an i nevitable conse- 
quen ce of the mass continuity (see Eq.9 in IBlinnikov et al.l 
l2012j . and its discussion). One can consider both shocks 
as "glued together" into dense shell. Moreover, all the ki- 
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netic energy of the matter accreted onto the dense shell 
is radiated away. That is why we call this structure cold 
dense shell (hereafter CDS). Nevertheless, the region of the 
cloud in front of the CDS will be heated up and shine in 
narrow spectral lines. Depending on t he optical depth of 
the heated cloud, C DS may shine out l|Chugai et al.ll2004l ; 
IWooslev et ail 120071 ). In this case the observed spectrum 
may have narrow lines on a broad pedestal. In addition, the 
matter of the cloud will be accumulated in the CDS, grad- 
ually breaking it, and eventually the neglect of the relative 
shock velocities becomes incorrect. In futu re evolution one 
has t o investigate the stability of CDS f see Ivan Marie et al.l 
l2010h . 

The calculations show that under certain conditions the 
photosphere is close to the CDS and moves together wit h 
it (e.g., 2006gy. IWooslev et all [20071 ; iBlinnikov et alj|2012j ). 
One may find the matter velocity from the broad spectra 
features which determine the photosphere velocity, v p h = 
dR ph /dt. 

After that a set of models is built, which well describes 
the time dependence of the flux, the color temperature, the 
cloud and ejecta velocities and the cloud density. Such mod- 
els do not necessarily correspond to the same distance to 
the star. DSM will "scale" the set of models to the desired 
distance by changing the radius of the photosphere. 

Let the initial radius (unknown to us) is Ro, and Ri = 
Ro + ARi for i = 1, 2, 3, . . ., where ARi is already known 
from dR integrated over time. We assume here that the ob- 
servations are sufficiently frequent to allow measurements 
of the increments in radius dR p h = v p hdt for a number of 
points, where dt is a difference of time of the successive ob- 
servations. 

Then we use Eq. 4 from (jBlinnikov et al]|2012l ): 



( vi (Ro + ARi)yJ-KB u {T cui ) = 10 a: 



(1) 



Here £„i is the dilution factor, T c „; is the color tem- 
perature, obtained from the spectrum, function B is the 
blackbody intensity, F V i is the observed flux, and A v is the 
extinction in stellar magnitudes. All quantities are defined 
for the frequency v. And finally, D is an unknown distance 
to the star. 

A good model gives us a set of the T c „; for all ob- 
servational points. From the measured F V i and ARi we can 



find Ro and the combination a„ 
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D by the least 



square method. Instead of v one may use index s labelling 
one of the broad-band filters. To find the distance D we need 
to know A 3 e.g. from the astronomical observations. 

As we have noticed in (|Blinnikov et al.ll2012l ). values 
of £ do not vary too much from model to model and do not 
strongly depend on the photospheric radius. In the follow- 
ing we assume the correction factor £ = 1, which is close 
to the values of £ with a ccuracy about ~ 10% found in our 
radia tion hydro models |Wooslev et al.l [20071 : iMoriva et al.l 
2013) for the growing part of the light curve. Of course, 
an accurate modelling requires building a hydrodynami- 
cal model not only for the light curve but also for spec- 
tral line profiles with an account of dilution and projection 
effe cts as it has been done for the recent Cepheid mod- 
els dGautschvlll987l; ISabbey et al.lll995l: IStorm et al.ll201ll ; 



3 DISTANCE TO THE SN 2009ip 

We have taken observ a tional data from ( Prieto et alj 20121: 
Mauerhan et al.l 120121 : iPastorello et al.l 120121 ) and iPrietol 
20121 ) web page. 

Supernova SN 2009ip began to brighten rapidly after 
September 24 possibly due to i nteraction of the ei ecta with 
the circu mstellar material (see Prieto et alj (|2012l ) and dis- 
cussion m iMauerhan et all l|2012l )). Unfortunately, starting 
from September 28 the broad lines mostly disappeared possi- 
bly due to the increasing opacity of the heated circumstellar 
matter, and it becomes difficult to determine the velocity of 
the CDS (hereafter "shell" = CDS). 

We have used a short 2-day period after September 24, 
when one can determine the velocity of the shell. In ad- 
dition, at this period the lumino sity is proportional to the 
square of time ijPrieto et al.ll2012h . which corresponds to the 
constant expansion velocity and color temperature of the 

photosphere. 

Photometric data for SN 2009ip in R band jBesscll 
2005) of all three aut hors agree well (see footnote 28 in 
Pastorello et al.1 J2012l)). We h ave taken the first 36 points 



els (ll iautscny iy»/'; Sabbcv et al.111^51: IStorm et al 
iRastorguev &: Dambisl 201lT " Rastorguev et al. 20121 ) 



of the tabltFI from lPrietol (|2012l ) web page. These points cor- 
respond to the days which are most interesting for us. 

To estimate the expansion velocity v of the shell we 
have used the data obtained from the H a absorption com- 
ponent. At the epoch of a sudden in crease of the lum i nosit 
of S N 2009ip (~ 23-24 Sep tember) IPastorello et al.1 (1201 
and IMauerhan et al.l (|2012l ) indicate values v « 13000 km/s 
and « 13800 km/s, respectively. Of course, a full analysis 
needs to address the question, which is not yet considered: 
are the layers of hydrogen excited to the third level (which 
produced H a ) close to the photosphere. We have taken the 
matter velocity equal to v ~ 13400 km/s. 

W e have adopted an extin ction of Ar — 0.051 mag fol- 
lowing [MauerhanTeTlil] ((2012) . In view of the smallness of 
this value, 10 2A R is close to unity. The error in Ar has 
very little effect on the final answer. We have taken 0.001 
mag for the error value and did not evaluate the impact of 
its change s on the final result. 

iPrieto et al.l (|2012T l obtains black-body temperature in 
two different ways. The first black-body fit takes into ac- 
count the data from R and I filters and has a value 
T ~ 14500 K. The second fit uses the longer wavelength 
coverage (fro m neax-UV to V-b and) and the temperature 
T « 19200 K. IPrieto et alj (|2012l ) treats the obtained values 
as effective temperatures, but the method of obtaining them 
corresponds to color temperatures. In our method these val- 
ues are used as color temperatures. 

To estimate the confidence intervals of the distance we 
have done a resampling Monte Carlo (MC) simulation based 
on these data. We resampled the values of T, of the stellar 
magnitude tur in standard filter R, the reddening Ar, and 
velocity v, each with normal distribution. For some of those 
quantities we have taken standard deviations from avail- 
able data or fixed "manually" for illustration if they are 
not known. For example, since we do not know the error of 
v, T, we examine how variations in their errors affect the 
response. 



1 http : //www. astro .princeton. edu/~ jprieto/sn09ip/photR. dat 
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Figure 1. Monte Carlo resampling simulation of the distance 
D to SN 2009ip by DSM with respect to temperature T. The 
isocontours of probability distribution functio n (pdf) are shown 
with equal step in pdf . The observations from jPrieto et al.l[20l3 : 
iMauerhan et alj|20lj; iPastorello et alj|2012ft have been used for 
36 different time points from |PrietcTi(2012h web page. 



For obtaining the confidence intervals for the mean and 
median distances it was sufficient to do 10 5 MC test s. 

Taking the temperature of the first IPrieto et all (|2012h 
estimate as T ~ 14500 K and the relative errors in v and 
T equal to 5% we have obtained both the mean and me- 
dian distances D ~ 16.1 Mpc with 68% confidence interval 
±0.6 Mpc. 

If we artificially increase by factor 2 the relative errors of 
v, T we find the distance D « 15.1 Mpc with 68% confidence 
interval (—1.1, +1.2) Mpc. The answer is more sensitive to 
the error in the temperature. 

T aking the temperature of the second IPrieto et al.l 
|2012ri estimate as T ~ 19200 K and the relative error 
in v and T equal to 5% we have obtained the distance 
D « 20.1 Mpc with 68% confidence interval ±0.8 Mpc. 
For illustration of this case, we have built a slice of multi- 
dimensional probability density functions in the plane - T 
and D. The plot in Fig. [1] is built with 10 7 samples to ob- 
tain a better statistics near the top of the distribution. The 
largest thick-line contour in Fig. [T] is about one standard 
deviation. 

If we increase by factor 2 the relative errors of v and 
T we find the distance D m 19.9 Mpc with 68% confidence 
interval (—1.4, ±1.5) Mpc. The answer is more sensitive to 
error in the temperature again. In all the cases the error 
growth leads to smaller median and mean distances. 

Those experiments show that the results are rather ro- 
bust given the level of accuracy of the data and models. 

The distance modulus for the galaxy NGC 725, 
where the SN 20 09ip is situated, i s generally accepted as 
= 31.55 mag l|Smith et al.l l201oh which corresponds to 
20.4 Mpc. More detailed and complete studies (like SEAM) 
without the black-body approximation but with an assis- 
tance of a complete hydrodynamical model of the whole 



sequence of SN 2009ip outbursts can accurately determine 
the temperature. However, already now we can say that 
in the frameworks of our models the temperature estimate 
T ~ 14500 K is unsatisfactory (if one believes in generally 
accepted value of the distance). The goal of the current pa- 
per is not to obtain the most reliable distance but to show 
that it can be easily obtained with reliable data. Moreover, 
increasing the number of the observations would reduce the 
error in the final result as the square root of the number of 
observations. 

Since the host galaxy NGC 725 is relatively close, the 
Hubble parameter based only on the redshift z can be 
wrong because of a large peculiar velocity of the galaxy. 
On the other hand, based on the generally accepted value 
Ho = 71 km/s/Mpc the resulting value of the distance 
to the galaxy NGC 725 can be significantly different from 
the distance obtaine d by other methods. For instance, 
iLevesaue et al.l (|2012T ) uses the distance to NGC 725 which 
is equal to 24 Mpc. 

Our distance value is obtained by the new direct method 
and does not rely on the Cosmic Distance Ladder! It is 
needed to investigate the role of variations of the correction 
factors in different SN 2009ip models to check the robust- 
ness of our results. We present here the values for D only 
for the illustration of the efficiency of the method. 

4 CONCLUSIONS 

We have obtained the distance to SN 2009ip equal D ~ 
20.1 ±0.8 (68% CL) Mpc. The rel ative error is much sma ller 
than in the case of SN 2006gy (|Blinnikov et all 120121 1. If 
the observations of a supernova are of good quality, and 
has a large number of data points, the internal error of our 
method for such supernova will be very small. It is needed 
to develop a complete supernova SEAM-like model for the 
comprehensive work of the method, but with complex hydro- 
dynamics taking into account the entire sequence of previous 
outbursts. One should evaluate the correction factors (e.g. 
dilution factor) to eliminate possible systematic errors. 

The distance we found was based on the spherical sym- 
metric model of the explosion. The fact that it coincides with 
the generally accepted value shows that the initial epoch of 
the growth of the SN 2009ip luminosity is well described by 
an expanding spherical shell. 

Constraining of the cosmological parameters and the 
understanding of Dark Energy depend strongly on accurate 
measurements of distances in the Universe. Our results on 
SN 2009ip supports the conclusion that SNe Iln may be used 
for cosmology as primary distance indicators with the new 
DSM method. 
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